Distribution of resonances in the quantum open baker map 
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We study relevant features of the spectrum of the quantum open baker map. The opening consists 
of a cut along the momentum p direction of the 2-torus phase space, modelling an open chaotic 
cavity. We study briefly the classical forward trapped set and analyze the corresponding quantum 
nonunitary evolution operator. The distribution of eigenvalues depends strongly on the location 
of the escape region with respect to the central discontinuity of this map. This introduces new 
ingredients to the association among the classical escape and quantum decay rates. Finally, we 
could verify that the validity of the fractal Weyl law holds in all cases. 

PACS numbers: 05.45.Mt, 03.65.Sq 



I. INTRODUCTION 

Recently, there has been an upsurge of the interest 
in open quantum systems, whose properties are still less 
known compared to those of the closed ones. Besides 
their general fundamental importance, they are also of 
the utmost relevance in very active fields. We can men- 
tion a few of them, such as the study of the quantum 
to classical correspondence [l|, quantum dots 01, mi- 
crolasers having chaotic resonant cavities [3, 0,111 and 
chaotic scattering [g, 7] . These systems are characterized 
by a nonunitary quantum evolution. The corresponding 
operators have a set of right and left decaying nonorthog- 
onal eigenfunctions associated to them. Their complex 
eigenvalues Zi (also referred to as resonances) fall inside 
the unit circle when represented in the complex plane, 
i.e. vf — — exp(— Ti) < 1. The exponent T^ > is 
the usually called decay rate. 

The classical phase space of open chaotic systems is 
characterized by fractal sets associated with trajectories 
that remain trapped for infinite times. Those orbits that 
stay forever in the future define what it is called the for- 
ward trapped set, and those that stay forever in the past 
define the backward one. An initial classical probabil- 
ity uniformly distributed in the phase space decays at an 
exponential rate. This allows to define the so-called clas- 
sical escape rate jci- The intersection of both sets, that 
is, the set of trajectories which do not escape to infinity 
either in the past or in the future is called the repeller. 

Regarding the eigenstates, we can distinguish between 
short-lived and long-lived ones. The former (F^ ^ 1) 
are associated with the trajectories that escape from 
the system before the Ehrenfest time, while the latter 
(Fi = 0(1)) are related to the classical trapped sets, thus 
carrying the most relevant classical information. One 
of the most important properties of open quantum sys- 
tems is the conjectured fractal Weyl law. This law re- 
lates the mean density of resonances, h and the structure 
of the classical phase space. It predicts that the num- 
ber of long-lived states goes as ^ ft"'''^^', where d 
is a fractal dimension of the classical strange repeller. 



This law has been checked for a three disk system [8| 
and some quantum maps d, [l^, EH, El, El, E3| , and it 
is still being tested. But much less is known about the 
distribution of the resonances. There are some results 
obtained for random matrix models [ll|. Also, a scal- 
ing property has been numerically verified for the open 
kicked rotator [1]. Essentially, the classical escape rates 
of this system determine the quantum decay rates associ- 
ated with the long-lived eigenstates. However, there are 
no analytical results for the semiclassical limiting distri- 
bution. We shed new light into this open problem by 
concentrating on the spectral behavior of the most sim- 
ple models of open chaotic dynamics, i.e. open piecewise 
linear maps. Discontinuities are an essential part of these 
systems, being responsible for their chaoticity. Then it is 
very interesting to study their influence on the spectral 
behavior. 

In this work we focus on the quantum open baker map, 
which is a chaotic transformation of the unit square (2- 
torus) phase space. This is a paradigmatic model in clas- 
sical and quantum chaos and also in statistical mechanics 
[1, El]- Its relevance both, in fundamental studies and in 
applications to a wide range of areas, is difficult to overes- 
timate. As such, a deep knowledge of its spectral features 
is very important. We have investigated the behavior of 
the distribution of its eigenvalues as a function of the lo- 
cation of the escape region in phase space. We have found 
that the central discontinuity of this map plays a crucial 
role in the behavior of the spectrum. The quantum ef- 
fects can be related to the classical behavior, which we 
study very briefiy in order to support our explanations. 
But there are also important features of purely quantum 
character. In fact, the link between the classical escape 
and the quantum decay rates is more subtle for openings 
that overlap with the central discontinuity of our map 
than for those which do not. The shortest periodic or- 
bits become relevant in the overlapping cases and based 
on this we provide a conjecture in order to explain our 
finding. 

This paper is structured as follows: in Section II we 
describe the model, giving a short introduction to the 
classical dynamics and the quantization method used. In 
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Section III we first study some aspects of the classical dy- 
namics that help us to understand the spectral behavior. 
Then we show the results for the distribution of eigenval- 
ues and present a conjecture explaining them. We also 
verify the validity of the fractal Wcyl law for all cases 
under investigation. Finally, in Section IV we draw the 
conclusions. 



II. THE CLASSICAL AND QUANTUM OPEN 
BAKER MAP 



In this Section we introduce the main features of two 
dimensional torus open maps and define our system. We 
use a very simple method to model dissipation as it oc- 
curs in scattering through a cavity, for instance. We 
assume that all the classical initial conditions that are 
mapped inside the area of phase space corresponding to 
the opening leave the system. Thus, the open map is 
defined on a subset of the 2-torus. If we choose an open- 
ing that represents a fraction M/N of the total area in 
phase space, then the open quantum map corresponds to 
a nonunitary matrix B" — BP. In this expression, P 
is a projector onto the complement of the opening, and 
the operator B corresponds to the closed quantum map, 
a unitary matrix acting on a Hilbert space of dimension 

N = l/2TTh. 

We study the open baker map with different opening 
sizes along the p direction. The closed classical trans- 
formation is defined in the 2-torus = [0, 1) x [0, 1) 

by 



with 



B{q,p) = 



(2g,p/2) 

(2(7-l,(p+l)/2) 



if < q < 1/2 
if 1/2 < g < 1 



(1) 



This transformation is an area-preserving, uniformly 
hyperbolic, piecewise-linear and invertible map with Lya- 
punov exponent A = In 2. Geometrically, the map 
stretches the unit square by a factor of two in the q di- 
rection, squeezes it by the same factor in the p direction, 
and then stacks the right half onto the left one. The 
opening is performed by eliminating from the evolution 
those initial conditions falling inside a rectangle of width 
Aq, centered at q^ and extending along the whole p axis 
in phase space. The dynamics of the open baker map 
has been previously studied [ll,[l^. However, the role 
played by discontinuities (specially the one along the line 
q = 1/2) has not received much attention. Moreover, 
we do not know of any study regarding its effects on the 
spectral behavior of the quantum version. 

Following the quantization process described in [l^, 
[iTj . in an even iV-dimensional Hilbert space, the quan- 
tum baker map is defined in terms of the discrete Fourier 
transform in the position representation as 



(Gjv)j7c 



exp{-27ri(j + l/2)(fc + l/2)/A^}. 



Bn — G 



N 



Gjv/2 

Gn/2 



(2) 



B^ is a unitary matrix and represents the quantum 
dynamics of the closed baker map. A Aq wide cut is 
made along the p direction by means of the projector op- 
erator P on its complement. Finally, the corresponding 
quantum dynamics for the open baker map is given by 
the nonunitary matrix = Bj^P. It is worth men- 
tioning that we have chosen to open the baker map in a 
different way than in [ill, [13 • This allows us to vary the 
location and width of the escape region. 



III. THE EIGENVALUES: THE ROLE OF 
DISCONTINUITIES 

A. The classical repeller 

We first study some features of the classical phase 
space of our system. This is solely intended to under- 
stand which are the main classical ingredients that play 
a significant part quantum mechanically. These ingredi- 
ents will help in the explanation of the quantum behavior, 
which is the main interest of this work. For that purpose, 
we have calculated the escape rate for different locations 
and sizes of the opening. By means of the area of the 
forward trapped sets as a function of the number of the 
iterations of the map Afw{t)^ the classical escape rate 
can be easily calculated as 7c/ — — \TiAfw{t)/t. Then, 
the information dimension di of the corresponding re- 
peller can be determined through the known relationship 
d/ = 2-7c,/A [6, 20, 21]. 

In Fig. [1] we can see the value of A/u,(10), i.e. the 
area of the tenth iteration of the map as a function of 
the position of the center of the opening qc, for three 
different values of the width, Aq = 0.05,0.1 and 0.2. 
Given that this quantity is symmetrical with respect to 
qc = 0.5 (for a given Aq), we only show values in the 
range qc G [0;0.5]. In order to calculate these curves 
we have evolved initial conditions uniformly covering the 
phase space, eliminating the area corresponding to the 
opening at each iteration. 

It can be clearly noticed that there is a common shape 
regardless of the size of the escape region. The minimum 
of is at qc = 0.5, while the maximum is generally 
at — 0. Also, we can identify a minimum at qc ^ 
0.25 and a maximum at qc ^ 0.3. This can be roughly 
explained by means of the first iterations of the opening 
through the map. In fact, the openings located at the 
central discontinuity typically do not overlap with their 
first iterations since there are no periodic orbits at g = 
0.5. As a result, initial conditions escape faster. This is 
not the case when the cut is made at the discontinuity 
at g = 0, which includes the shortest periodic orbit (q = 
0,p = 0). As a consequence, an opening overlapping 
with this discontinuity behaves much in the same way 
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as one having a generic qc value. We have found that it 
is possible to use the shortest periodic orbits to give a 
good estimate for the main features of these curves (this 
will be explained elsewhere [13] )■ In the insets we can 
see two examples of the shape of the forward trapped 
sets at qc = 0.3 and qc = 0.5 for Aq = 0.05. These 
two positions of the escape region illustrate the two most 
relevant situations: in the first case the opening is far 
from the central discontinuity of the map and its first 
iterations, while in the second case it is centered on it. 
The maximum escape rate is reached for this latter and 
this can be easily associated with the thinner look of the 
phase space distribution. We have calculated the escape 
rates for Aq = 0.1 with qc = 0.3 and qc = 0.5, resulting 
in 7c/ = 0.09073 and 7^; = 0.16488, respectively. The 
corresponding information dimensions are di = 1.86910 
and dj — 1.76213. In the following we will see how this 
behavior translates into the quantum domain. 




'0 0.1 0.2 q 0.3 0.4 0.5 



Figure 1: (Color online): Area of the forward trapped set for the 
10th iteration of the open baker map Af^{10) as a function of the 
center of the opening qc for different widths, with Ag = 0.2 (black 
solid line), Ag = 0.1 (red dotted line) and Aij = 0.05 (green dashed 
line). The insets show the forward trapped sets (in black) for t = 10 
and A<j = 0.05, for <jc = 0.3 (left) and = 0.5 (right). 



B. The eigenvalues 

We now study the behavior of the distribution of the 
eigenvalues of the quantum evolution operator. We first 
show the eigenvalues in the complex plane for Ag = 0.1. 
As can be seen in Fig. [2l since moduli are less than one 
[v <1) all of them fall inside the unit circle. They clus- 
ter near the origin and at a ring close to = 1. In the 
upper panels N — 602, while in the lower ones N — 2048. 
In both cases we display the results for qc — 0.3 on the 
left and qc — 0.5 on the right. In these last situations 
there is a much less dense distribution of eigenvalues at 
the outer ring, and an increase of density near the origin. 
This is consistent with the predictions of the fractal Weyl 
law, given that the information dimension of the classical 
repeller is smaller for qc — 0.5. It is interesting to men- 
tion that for an opening at gc = the behavior is rather 



similar to what happens for one at a qc far from discon- 
tinuities. To gain further insight about the behavior of 




Figure 2: Eigenvalues of the open baker map in the complex 
plane. In the upper panels we show the case for A'^ = 602 and 
in the lower ones for A'^ = 2048. In the left panels the opening 
is centered at qc = 0.3 and in the right ones at qc — 0.5; 
Aq = 0.1 in all cases. 



this distribution, we have calculated the normalized cu- 
mulative number of resonances n = i /N as a function of 
V. In Fig. [3] results for the same values of Ag, TV and 
qc than in Fig. [2] can be seen (see caption for details). 
From this Figure it is clear that there is a higher den- 
sity of eigenvalues near v = Q for qc = 0.5, and a lower 
density near v — 1 for qc ~ 0.3. Also, the shape of the 
tails [v > 0.7) of both distributions seems to be different, 
showing a more linear behavior in the former rather than 
in the latter case. 

In order to better evaluate the eigenvalue distribu- 
tion we have constructed the histograms for W = dn/du 
where the bin size has been taken /S.v = 0.01. We show 
values for v > 0.7, corresponding to the tails of Fig. [3] To 
compare the distributions W at different relevant cases, 
we have constructed Fig. 21 In these plots we have su- 
perimposed the cases for qc — 0.3 and gc = 0.5. In the 
upper panel the N — 2048 case can be seen. This is an 
example of the special situation for the baker map when 
the dimension is of the form N = 2^ where I is an inte- 
ger number. In the middle panel we can see an example 
for N = 602 (in this case N ^ 2^) where the width for 
qc — 0.5 is much greater than that for qc = 0.3. In the 
lower panel there is an example for N — 1782 that shows 
a smaller width, now for the qc — 0.5 case. 

It is clear from Fig. 2] that the width of W can vary 
significantly with TV. To have a complete picture of 
this behavior we have also computed the width a of the 
eigenvalue distributions as a function of the dimension 
of the Hilbert space in the interval N e [500; 2000], for 
Aq — 0.1. We have numerically measured the width of 
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Figure 3: (Color online) Cumulative number of resonances n 
as a function of v. The green short dashed line corresponds 
to iV = 2048 and the black dot dashed line to iV = 602, both 
for Qc = 0.5; the blue long dashed line to N = 2048 and the 
red dotted line to A'^ = 602, both for qc = 0.3. 
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Figure 5: (Color online) Width g of the eigenvalue distribu- 
tions W as a function of the dimension of the Hilbert space. 
The red dotted line corresponds to qc — 0.5, and the black 
solid line to qc = 0.3, being Aq — 0.1 in both cases. Only 
values for v > Q.7 were considered. 




Figure 4: (Color online) Histograms corresponding to the eigen- 
value distribution W as a function of v for qc = 0.3 (red dot- 
ted lines) and Qc = 0.5 (black dashed lines). In the upper panel 
A'^ = 2048, in the middle panel A'^ = 602, and in the lower panel 
N = 1782. 



each liistogram at half height for v > 0.7. Results are 
shown in Fig. [Sj where we have taken qc = 0.3 and 0.5. 
We can see that the widths of the distributions are sim- 
ilar for a wide range of Hilbert space dimensions. They 
are generally greater for qc — 0.5 than for qc — 0.3 by 
approximately a factor of 1.5. However, there are peaks 
for specific values of N where the opening contains the 
central discontinuity, that show both much higher and 
lower widths than for the qc = 0.3 case. In fact, we can 
identify the two peaks corresponding to the cases shown 
in the middle and lower panels of Fig. [H^i.e., N — 602, 
and N = 1782). It is interesting to note that these peaks 
seem to be present also at the high N limit. 

To see more clearly the scaling property of the distri- 
butions shown in Fig. [51 we have represented the same 



data, but now as a function of the decay rate and rescaled 
with the classical escape rates 7c/. We show this in Fig. 
[S] From it, we can see that for openings with qc = 0.3 
all distributions are almost the same, with a clear peak 
falling at r ~ 0.1 ~ 7c/. On the other hand, openings 
overlapping with the central discontinuity show eigenval- 
ues distributions that are not rescalable to the ones cor- 
responding to the previous non-overlapping cases. More- 
over, they are not rescalable among themselves. Fluc- 
tuations become the rule in these cases, ranging from 
peaks narrower than in the generic situations (like for 
N — 1782) to distributions where no clear maximum can 
be found (hke for N = 2048 and N = 602). This comple- 
ments previous results found in the literature Q for the 
open kicked rotator, where this scaling turned out to be 
universally valid. 

We could associate this behavior with the fact that 
in the qc — 0.5 case the shortest periodic orbits survive 
the dissipation process and have a substantial role in the 
localization properties of resonances. In fact, there is 
numerical evidence that individual resonance eigenstates 
of an open quantum system present localization around 
unstable short periodic orbits in a similar way as their 
closed counterparts [14]. This so-called scarring phenom- 
ena could be important enough to make almost dissapear 
the classical escape rate information from the quantum 
distributions. This last value is an average and now fluc- 
tuations become relevant. This includes situations where 
some particular quantization condition shrinks the distri- 
bution for a given N, as in the case of iV = 1782. In this 
new situation, the single quantum decay rates of reso- 
nances associated to given short periodic orbits could be 
singled out, and this would explain the multipeak struc- 
ture of the corresponding distributions. When many or- 
bits are relevant, and also longer ones, their decay rates 
combined could be more easily associated to a value that 
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Figure 6: (Color online) Same distributions as in Fig. |4] shown as 
a function of the decay rate and rescaled with 7c! ■ Thin red lines 
correspond to Qc = 0.3 and thick black ones to qc = 0.5. Solid lines 
correspond to N = 2048, dashed to N = 602 and dot-dashed to 
N = 1782. 



should better approximate the classical one. 

Finally, we have checked whether the fractal Weyl law 
is verified or not by three different widths of the opening 
and for both representative Qc values. The logarithmic 
plots of the fraction of eigenstates for > 0.3 as a func- 
tion of the dimension of the Hilbert space can be seen 
in Fig. [71 The lines correspond to the prediction of the 
fractal Weyl law log (N^) = log {N){di - I) + A (where 
A is a constant). The symbols correspond to the numer- 
ically calculated values (see caption for more details). In 
all of these cases the slope is correctly described by the 
information dimension calculated from the area of the 
forward trapped set. We have adjusted the constant in 
order to fit the data. While in the left panel the central 
discontinuity is not inside the opening and in the right 
one it is, we find that in both cases the agreement with 
the theoretical prediction is very good. 



IV. CONCLUSIONS 

In this work we have investigated the behavior of open 
piecewise linear maps. We have found that in the open 
baker map the role played by the central discontinuity is 
crucial. Firstly, we have made a brief study of the proper- 
ties of the classical forward trapped sets. We have calcu- 
lated the escape rates as a function of the location of the 
opening with respect to this singularity. We have found 
that when a cut along the p direction contains the dis- 
continuity the information dimension of the repeller goes 
to the minimum. On the other hand, we have studied the 
behavior of the quantum map. We could verify that the 
distributions of the eigenvalues with the greatest classi- 
cal information (i.e., with the greatest v) show a similar 
behavior for different A'' values when the opening is lo- 




Figure 7: (Color online) Logarithmic plot of the fraction of eigen- 
states N,j for u > 0.3, as a function of the dimension of the Hilbert 
space N. Lines correspond to the prediction given by the fractal 
Weyl law. Symbols correspond to the numerically calculated val- 
ues. In the upper panel qc = 0.3 and in the lower panel qc = 0.5. 
In both cases black solid lines and circles correspond to Ag = 0.05, 
red dotted lines and squares to Aq = 0.1, and green dashed lines 
and triangles to Aq = 0.2. 



cated at a typical value of q (far from the influence of the 
central discontinuity). But when this is not the case the 
eigenvalues behave in a non standard way, showing dis- 
tributions that cannot be rescaled to the previous ones, 
and not even among themselves for different N. For the 
less dense fractals corresponding to a smaller dimension 
(Qc = 0.5) the role played by the shortest periodic or- 
bits becomes relevant. This supports our conjecture that 
they are responsible for this behavior. Also, fluctuations 
connected with the quantization rules for these trajecto- 
ries seem to be responsible for the great differences found 
for different values of N. Then, the results presented in 
this work could be of much relevance in order to under- 
stand the fundamental problem of scarring phenomena 
present in open quantum systems. We think that after 
developing a suitable probe of localization for this case, 
an interesting study will consist of checking if the peaks 
found in the W distribution correspond to short periodic 
orbits. Moreover, it will be very important to see if this 
behavior is also present in other types of systems, like 
open kicked maps, for example. We will investigate on 
this in future studies 22]. Finally, it is remarkable that 
the fractal Weyl law is obeyed with great accuracy in all 
cases. 
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